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Group velocity of 17.8-kc/s VLF radio waves is determined by measuring the difference in time 
of arrival of the transmitted signals which travel over complementary great circle paths. The velocity 
is found to be 3.04 X 10 8 meters/second with an estimated error of approximately ±0.10 X 10 8 meters/ 
second. Differences in measurement technique and general environmental conditions between this 
experiment and a previous measurement of 18-kc/s group velocity which was made in 1949 is discussed. 



1. Introduction 

1.1. Equipment and Timing 

In a previous article [Westfall, 1964], a system for 
simultaneously measuring phase and amplitude of 
VLF radio waves over complementary great-circle 
paths was described. Sample data for the long and 
short great-circle paths were shown where the trans- 
mitter was NAA and the receiving site was at the U.S. 
Navy Electronics Laboratory at San Diego, Calif. 

The long-path receiver "looks" for the long-path 
signal during the off periods of the short-path signal. 
The amount of time delay between the arrival of the 
short-path signal and the arrival of the long-path signal 
at San Diego may be used to measure the group ve- 
locity of the VLF signals. While this parameter can- 



not be determined pulse by pulse because the long-path 
signal is quite weak, a coherent detection system pro- 
vides an integrated long-path signal voltage which will 
be a maximum when the delay between the arrival 
of the short-path signal and the "look time" on the long- 
path receiver is optimum. 

Figure 1 is a block diagram of the long-path portion 
of the receiving system which was used to make the 
group velocity measurements. 

The antenna system is composed of commercially 
available loop and whip antennas and appropriate 
circuitry to produce a cardioid receiving pattern with 
a null of about 30 dB. This null is normally directed 
toward the transmitter along the short great-circle 
path. The VLF phase-tracking receiver is also of 
commercial design and provides an analog record 
of VLF phase. The phase coherent 1 kc/s IF is passed 



20 6 
FOOT FOOT 
WHIP LOOP 
ST7 J=Z) 














RUBIDIUM 
FREQUENCY 
STANDARD 








SERVO CONTROL 


d-c 
AMPLIFIER 




\ 


' Li 


JJ 


























w 
















CARDIOID 
UNIT 








VLF PHASE 

TRACK 

RECEIVER 




DELAY 
LINEI 




CONTROL 
GATE 




PHASE 
DETECTOR 










IF 




































L_ 


REFERENCE 

kc/s <o° 


10 
FOOT 
WHIP 


LONG PATH 
PHASE 










PHASE 
DETECTOR 




LONG PATH 
AMPLITUDE 


\ 


/ 












L_ 


REFERE 

kc/s 


NCE 

<90° 






AUXILIARY 

VLF 
RECEIVER 




GATE 

TRIGGER 

GENERATOR 






LONG 
PATH 






































ANALOG 




























RECORDS 

OF 
RECEIVING 




DELAY 
LINE 2 






DELAY 
LINE 3 






SHORT 
PATH 




TIME 

















































FIGURE I. VLV receiving system designed to reject the short-path signal and measure phase and 
amplitude of the long-path signal. 
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through delay line 1 and enters the gating circuits. 
The gating circuits are controlled by the VLF signal 
received via the short path on the auxiliary receiver. 
The gates are normally turned on only when the short 
path signal is off and also only at the approximate 
time of arrival of the long-path signal. This timing 
arrangement is accomplished by conventional pulse 
circuitry in the gate trigger-generator. 

Figure 2 shows this timing in somewhat idealized 
detail. Because of the slope in the leading and trailing 
edges of the short-path signal, Pulse A is narrower 
than the signal which produced it. The trailing edge 
of Pulse B is delayed by delay line 2 for a time interval 
marked T2. The short-path signal is delayed by delay 
line 1 for a time interval marked T\ until it coincides 
with Pulse B. During Pulse B, the receiver is turned 
off by the control gate. Pulse A is also delayed by 
delay line 3 until it coincides with the arrival of the 
long-path signal. At this time, the control gate turns 
the receiver on if no short-path signal is present. In 
this example, the short-path signal is so long that the 
long-path receiving gate is shorter than the long-path 
signal and only a portion of the long-path signal is 
received. 

The gate trigger generator also provides analog 
voltages which give the percent of time that the 
short-path signal is present and also the percent of 
time that the long-path receiver looks for the long- 
path signal. 

The integrated long-path amplitude voltage must be 
normalized because the look times are at random 
periods and the amount of time that the long-path 
receiver is actually receiving long-path signal will 
vary depending upon the nature of the code trans- 
missions. The percentage of time that the long-path 
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FIGURE 2. Long-path VLF receiver timing diagram. 



receiver is turned on is recorded as an analog voltage, 
and the normalized voltage V n is defined as 

v= v - 

Vn D 

where V is the integrated long-path signal voltage 
and D is the percentage of time that the long-path 
receiver is turned on. 

The delay time between short-path signal arrival 
and long-path "look" can be adjusted from zero to 
nearly 200 msec and an optimum arrival time for the 
long-path signal may be determined. 

2. Uncertainties Contributing to Measure- 
ment Error 

While the experiment appears to be simple, several 
factors combine to make the determination of a pre- 
cise value for group velocity quite difficult, i.e., 

(1) The random nature of code transmissions makes 
it difficult to determine precisely the average shape of 
the short-path keyed characters. It is necessary to 
know this shape in order to define a zero time. 

(2) The long-path signal amplitude may be constantly 
fluctuating in an unpredictable manner, which re- 
duces the significance of any single set of consecutive 
readings. 

(3) The amount of atmospheric noise, short-path 
signal amplitude, and subtle changes of wave shape 
in the short-path signal may influence the operation 
of the timing gate generator used in this system. 

(4) The delays inherent in the long-path receiver 
and the control receiver which provides signal to op- 
erate the timing gates are not the same and may not 
be constant under all conditions. 

(5) Finally, the long-path receiving gate is some- 
what shorter than the short-path signal which produces 
it. This has the effect of broadening the period of 
optimum time delay. 

3. Experimental Procedure and Results 

Figure 3 shows an estimate of the smoothed response 
curve throughout a 200 msec delay range. The data 
points were obtained by delaying pulse A (fig. 2) with 
respect to the time of arrival of the short-path signal. 
This figure is the result of about 400 individual read- 
ings taken during different hours of the day on 17 
days. Since the attenuation rate changes with time, 
the integrated amplitude voltage of the long-path signal 
varies with time, and as a result V n varies widely for 
the different sets of data. The value of V n which 
was nearest to 104 msec was made unity for each set 
and all values in the set were normalized accordingly. 
The response to small time delays is not shown because 
in this area D is relatively small and small errors in 
D produce large errors in V n . At zero delay, V equals 
zero. The data points in figure 3 suggest an optimum 
time delay, r opt , of about 105 msec with an estimated 
uncertainty of ± 5 msec which is due to the previously 
mentioned sources of experimental error. 
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A calibration of the actual delay, which is equal to 
the quantity T3 — T\ (fig. 2), was obtained by reversing 
the antenna cardioid pattern and photographing the 
short-path signal and the delayed pulse A while they 
were simultaneously displayed on a dual-beam oscil- 
loscope. To determine the actual delay time for a 
given delay line setting, the delays between the lead- 
ing edges as well as the trailing edges were measured. 
One half the sum of the averages of these delays rep- 
resents the average delay r avg for a given delay line 
setting. Figure 4 shows the measured r avg for delay 
line settings of 110, 120, and 126. 

Further knowledge of the possible uncertainty of 
this curve may be obtained if it is assumed that the 
deviations in time of the leading and trailing edges from 
their respective averages is independent from the 
actual delay line settings, and also that the standard 
deviation of the leading and trailing edges is approxi- 



mately equal to 



2(Ar) 2 



for either case, where At 



is the deviation from the average and n is the number 
of samples. For approximately 50 samples of each, 
the leading edge standard deviation was about 1.58 
msec, and the trailing edge standard deviation was 
1.55 msec. The combined standard deviation, S, 
which specifies the standard deviation for the overall 
delay is therefore 2.2 msec. From these considera- 
tions alone it appears that the line in figure 4 should 
have an upper and lower uncertainty of about 2 msec. 

A second and perhaps more accurate method of 
determining r opt is to open the servo control circuit 
shown in figure 1 and to apply the quadrature d-c volt- 
age outputs of the phase detectors to the vertical and 
horizontal inputs of an oscilloscope. The resulting 
deflection of the beam spot (which is a direct indica- 
tion of amplitude voltage V) may be recorded with 
suitable time exposure photographs. Repeated meas- 
urements with different delay line settings around the 
estimated optimum time of arrival determine a range 
of optimum delay setting values which cannot easily 
be reduced. 

On 26 February 1965 NAA departed from their usual 
FSK schedule and maintained unbroken code trans- 
missions for several hours. This permitted a series 
of consecutive photographs to be made of both the 
circuit delay as well as the optimum time delay T opt , 
thus improving the accuracy of this series of measure- 
ments. About 120 photographs were made in a period 
of 3 hr. The region in figure 4 marked "optimum 
range" is the result of the measurement technique 
just described. It is an area in which a maximum 
amplitude V n appeared to exist, but within which the 
point of peaking could not be exactly determined. 
From figure 4 it is concluded that T opt was 102 ±3.5 
msec. 

If it is assumed that the group velocity over the first 
4.5 megameters (Mm) of the long path is the same as 
that over the short path (which is 4.5 Mm long), the 
optimum time delay measurement applies to a path 
31.0 Mm long and the resulting computed group 
velocity is 3.04 X 10 8 m/sec with an estimated error 
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FIGURE 3. Average relative amplitude as a function of time delay. 
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FIGURE 4. Measured actual delay as a function of delay line setting 
and measured optimum range for maximum V n . 

of ±0.10X10 8 m/sec. If for physical reasons group 
velocities greater than the speed of light are rejected, 
the final limits applying to this measurement would 
be 2.94 X 10 8 m/sec to 3.0 X 10 8 m/sec. 

From published theoretical curves for phase vel- 
ocity [Wait and Spies, 1964] it appears that there is 
little difference between east-west and west-east 
daytime propagation at 17.8 kc/s and therefore the 
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group velocity would not be expected to differ ap- 
preciably for the two directions. Also, experimental 
data at 10.2 kc/s [Pierce and Nath, 1961] shows that 
the daytime phase velocities over land and sea differ 
by only about 0.05 percent which suggests that the 
travel time of the radio waves would not differ appre- 
ciably for the two 4.5-Mm paths. 

If a group velocity V g is expressed as 



v,- 



n 



W- 



dV B 



df 



where / is frequency and V p is phase velocity, the 
phase velocity curves of Wait and Spies [1964] may be 
used to obtain a theoretical group velocity at this 
frequency of about 2.98 X 10 8 m/sec. It should be 
noted, however, that both V p as well as the derivative 
of phase velocity with respect to frequency can change 
slightly with direction of propagation, ionospheric 
height, and ground conductivity. 

This result agrees with the above measured value 
within the experimental error. 



Diego except that the path over the continental United 
States added another 4.5 Mm of land to the 18-kc/s 
path. The amount of land encountered on the long 
path divided by the total path length for the 18 and 
17.8 kc/s measurements produce ratios of about 0.29 
and 0.21, respectively. 

It should also be observed that the 1949 measure- 
ments were made during the peak of the 18th solar 
cycle when the smoothed sunspot number was about 
130. The present solar activity is at a minimum. A 
recent article [Poppoff et al., 1964] points out that the 
Z)-region structure is controlled by x radiation, cosmic 
rays, and Lyman-alpha, all of which change throughout 
the solar cycle. 

The 18-kc/s measurements were made on 25 January 
and 8 March 1949 with a difference of approximately 
0.8 msec in delay being attributed to disturbed con- 
ditions on 25 January. The NAA group velocity data 
were obtained in October 1964 and February 1965. 
Seasonal effects are being observed in the NAA long 
path phase and amplitude records, but the manner in 
which these changes are related to possible group 
velocity changes must await further study. 



4. Previous Group Velocity Measurement 

An early measurement of the group velocity of 
special 18-kc/s pulse transmissions was made by Brown 
[1949]. Round-the-world signals were observed at a 
site close to the transmitter at Annapolis, Md. A 
large loop was directed east-west and on days when 
the atmospheric noise was low, the long-path signals 
were directly visible on an oscilloscope and could be 
photographed. The group velocity determined in this 
way (assuming a 40 Mm distance) is about (2.91 ± 0.01) 
X 10 8 m/sec. 

It is interesting to compare the conditions existing 
in 1949 with the conditions currently existing. It 
seems likely that the signals which Brown measured 
arrived primarily from the west because of antenna 
orientation and because attenuation rates for east-west 
and north-south propagation are greater. The NAA 
signal arrives from the west also. The 18-kc/s signals 
probably traveled over nearly the same path as the 
current NAA transmissions being measured at San 



The author is grateful to A. D. Watt and J. E. Bickel 
for helpful suggestions during the preparation of this 
material. 
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